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We use the so-called Liouville-von Neumann (LvN) approach to study the 
nonequilibrium quantum dynamics of time-dependent second order phase 
transitions. The LvN approach is a canonical method that unifies the func- 
tional Schrodinger equation for the quantum evolution of pure states and the 
D \ LvN equation for the quantum description of mixed states of either equi- 

librium or nonequilibrium. As nonequilibrium quantum mechanical systems 
we study a time-dependent harmonic and an anharmonic oscillator and find 
^ ■ the exact Fock space and density operator for the harmonic oscillator and 

the nonperturbative Gaussian Fock space and density operator for the anhar- 
monic oscillator. The density matrix and the coherent, thermal and coherent- 
thermal states are found in terms of their classical solutions, for which the 
effective Hamiltonians and equations of motion are derived. The LvN ap- 
proach is further extended to quantum fields undergoing time-dependent sec- 
ond order phase transitions. We study an exactly solvable model with a 
finite smooth quench and find the two-point correlation functions. Due to the 
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spinodal instability of long wavelength modes the two-point correlation func- 
tions lead to the i 1//4 -scaling relation during the quench and the Cahn- Allen 
scaling relation i 1 / 2 after the completion of quench. Further, after the finite 
quench the domain formation shows a time-lag behavior at the cubic power 
of quench period. Finally we study the time-dependent phase transition of a 
self- interacting scalar field. 
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I. INTRODUCTION 



A system can interact directly with an environment to make its coupling parameters 
depend explicitly on time. Even the effective coupling parameters of a subsystem of a closed 
system, though conserved as a whole, may depend implicitly on time through an interaction 
with the rest of the system. The characteristic feature of these open systems is that their 
effective coupling parameters depend explicitly or implicitly on time. Therefore the genuine 
understanding of these systems requires the real-time processes from their initial conditions. 
Of a particular interest are the systems undergoing phase transitions. When a system cools 
down through an interaction with an environment, it may undergo a phase transition and 
its coupling parameters depend explicitly on time. Similarly, matter fields in the expanding 
early Universe undergo phase transitions through the interaction with gravity. 

Phase transitions are one of the most physically important phenomena in nature and have 
wide applications from condensed matter physics, particle physics and even to cosmology. 
The Kibble mechanism explains formation of topological defects in symmetry breaking phase 
transitions |I| . The kinetic process by Zurek has revealed a new feature of symmetry breaking 
phase transitions 0. The dynamics of phase transitions and formation of topological defects 
since then have become an important tool in variety of phenomena in condensed matter 
physics ||. It is also widely accepted that symmetry breaking phase transitions in the 
early Universe are inevitable for the structure formation of the present Universe |4j]. There 
has been an attempt through laboratory experiments to investigate the process of structure 
formation in the early stage of Universe ||. In QCD the quark-antiquark condensate breaks 
chiral symmetry when the temperature of quark-gluon plasma is lowered ||f7| . 

However, the most difficult and subtle facet of phase transitions is to understand its 
dynamics and the formation process of topological defects. As emphasized above, systems 
do become out of equilibrium in general during phase transitions because their coupling 
parameters depend explicitly on time through the interaction with an environment (heat 
bath). Hence the nonequilibrium dynamics of phase transitions should differ significantly 
from the equilibrium dynamics. In finite temperature field theory one obtains the effective 
action for the system in a thermal equilibrium or quasi-equilibrium by calculating its quan- 
tum fluctuations about a vacuum || . But as the phase transition proceeds, the fluctuations 
grow and the stability is lost. The effective action, when extrapolated to the phase transition 
regime in a literal sense, has a complex value, the imaginary part of which is related with 
the decay rate of the false vacuum |{J. In this sense finite temperature field theory can not 
be directly applied to study symmetry breaking phase transitions. 

Schwinger and Keldysh introduced the closed time-path integral to treat properly the 
quantum evolution of quantum fields out of equilibrium from their initial thermal equilib- 
rium |10| . Since then the closed time-path integral method has been developed and applied 



to many related problems |JL1|| . Recently the closed time-path integral method has been 



employed to study the nonequilibrium dynamics of second order phase transitions [12 1 G 



Another method is the functional Schrodinger-picture approach, in which the evolution of 
quantum states is found for explicitly time-dependent Hamiltonians ||17||. The large N- 



expansion method fll8||l9l and the mean-field or Hartree-Fock method |20|-p2| are used in 
conjunction with either the Schwinger-Keldish or functional Schrodinger method. Still an- 
other methods are the time-dependent variational principle J23| , the generating function for 
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correlation functions | 24| , and thermal field theory f25 



The purpose of this paper is two-fold. In the first part of this paper we elaborate 
and establish the recently introduced Liouville-von Neumann (LvN) approach so that it can 
readily be applied to nonequilibrium dynamics. In the second part we apply the approach to 
the systems undergoing time-dependent second order phase transitions. The LvN approach 
is a canonical method that unifies the functional Schrodinger equation for the quantum 
evolution of pure states and the LvN equation for the quantum description of mixed states 
of either equilibrium or nonequilibrium. One of the advantages is that one can make use of 
the well-known techniques of quantum mechanics and quantum many-particle systems. It 



is based on the observation by Lewis and Riesenfeld |26j that the quantum LvN equation, 
which is originally used to define the density operator for a mixed state, can also be used to 
find the exact pure states of a time-dependent harmonic oscillator. This observation makes it 
possible to find not only the mixed state but also the pure state of a time-dependent system. 
This LvN approach has been developed and applied to quantum fields in an expanding 



Friedmann-Robertson- Walker universe [p7j— p9(1 and to open boson and fermion systems |30 

In this paper we particularly focus on the model systems whose coupling parameters 
change signs during the evolution and emphasize the role these systems playing in the second 
order phase transition. In the case of a time-dependent harmonic oscillator or an ensemble 
of such oscillators with a positive time-dependent frequency squared, the Fock space con- 
sists of the number states which are the exact quantum states of the Schrodinger equation 
| 26| , |3T|| . The density operator is constructed in terms of the annihilation and creation oper- 



ators that satisfy the LvN equation pq| . Hence the coherent, thermal and coherent-thermal 
states are found rather straightforwardly according to the standard technique of quantum 
mechanics. We further show that the same construction of the Fock space and density op- 
erator still holds for the time-dependent harmonic oscillator with a sign changing frequency 
squared. By studying some analytically solvable models we investigate how the instability 
grows. Another technical strong point of the LvN approach is that it can also be used for a 
time-dependent anharmonic oscillator. At the leading order the LvN approach is equipped 
with the time-dependent annihilation and creation operators, the vacuum state of which is 



already the time- dependent Gaussian state that minimizes the Dirac action [33,^]. The 
LvN approach thus provides one with a nonperturbative quantum description for the time- 
dependent anharmonic oscillator, too. We find the coherent, thermal and coherent-thermal 
states for the anharmonic oscillator with a polynomial potential and study the dynamics 
from their effective actions. 

As field models for the second order phase transition, we consider first a free massive and 
then a self-interacting scalar field, the mass of which changes the sign in a finite time period 
through an external interaction. By studying analytically the exactly solvable model of the 
free scalar field, we show how the instability of long wavelength modes grows in time. The 
two-point thermal correlation function is expressed in terms of the classical solution for each 
mode that is already found in terms of a well-known function. The domain sizes are evaluated 
analytically by the steepest descent method and are shown to grow as t 1 / 4 during the quench 
and as the Cahn- Allen scaling relation t 1//2 after the completion of quench. Remarkably, the 
Cahn- Allen scaling relation shows a time-lag given by the cubic power of the quench period, 
which is absent in the instantaneous quench model [|12| , |r3f . However, both the instantaneous 
and the finite quench models have the same scaling relation given by the classical Cahn- Allen 
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equation ||, confirming the result of the instantaneous quench model |12|,13|. The free scalar 



field model describes only the stage of spinodal instability from the unstable local maximum 
to the spinodal line. In the self-interacting field model, the back-reaction of an interacting 
term contributes positively to the frequencies of long wavelength modes and shuts off the 
exponential growth of instability after crossing the spinodal line at sufficiently later times. 

The organization of this paper is as follows. It mainly consists of two parts: in the 
first part from Sec. II to Sec. V the LvN approach is elaborated to be applicable to phase 
transitions, and in the second part of Sees. VI and VII it is then applied to the second 
order phase transitions. In Sec. II we review the LvN approach to time-dependent quantum 
systems. In Sec. Ill we study the time-dependent harmonic oscillator and find the exact 
Fock space and the density operator. The density matrix is found and the nonequilibrium 
quantum dynamics is studied for the coherent, thermal and coherent-thermal states. In 
Sec. IV the LvN approach is applied to time-dependent inverted oscillators as a toy model 
for the second order phase transition. In Sec. V we extend the LvN approach to time- 
dependent anharmonic oscillators with polynomial potentials. The LvN approach leads to 
the nonperturbative Gaussian state at the leading order. We also show that the coherent 
state of the LvN approach recovers the result from the time-dependent mean-field or Hartree- 
Fock method. In Sec. VI we study an exactly solvable model of a free scalar field, which 
has a time-dependent mass coupling parameter and undergoes smoothly the second order 
phase transition for a finite quench period. The two-point thermal correlation function is 
analytically evaluated during the quench and after the completion of quench, and the scaling 
relations for the domain size are found. In Sec. VII we study a self-interacting scalar field 
that undergoes the time-dependent second order phase transition. 



II. LIOUVILLE-VON NEUMANN (LVN) APPROACH 

In this section we briefly review but emphasize the underlying assumptions of the LvN 



approach to time-dependent quantum systems introduced in Ref. [^7]. A time- dependent 
system can not remain in the initial equilibrium or the instantaneous quasi-equilibrium 
because the density operator 

,„ = ^f™\ (2.1) 

where H{t) is a time-dependent Hamiltonian operator and is the partition function, 
does not satisfy the quantum LvN equation. This means that though the system starts in 
the initial thermal equilibrium, its final state can be far away from the initial one. Even 
an initial pure state evolves toward a final one which differs drastically from the initial 
one, and leads, for instance, to particle production |33|-p^1. The Fock or Hilbert space 
of such time-dependent system of a quantum field transforms unitarily inequivalently so 
that even the initial vacuum state evolves to a superposition of particle states at final 
times. Thus the nonequilibrium system described by time- dependent quantum Hamiltonian 
follows the evolution of a mixed state that is out of equilibrium and is characterized by time- 
dependent processes. To properly describe the nonequilibrium evolution we shall adopt two 
assumptions. 
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First, from a microscopic point of view we assume that the nonequilibrium system obeys 
the quantum law, i.e. the time-dependent Schrodinger or Tomonaga-Schwinger equation 

zh^(t)) = H(t)Mt))- (2-2) 

Here H(t) is the time-dependent Hamiltonian of the system. This assumption is physically 
well-grounded since all the individual constituents of the system should obey the quantum 
law and the system as a whole should still obey the quantum law provided that all the 
interactions among the individuals are properly taken into account. 

Second, from a statistical point of view it is assumed that even the nonequilibrium system 
obeys the time-dependent quantum Liouville-von Neumann (LvN) equation 

ih^p{t) + \p(t),H{t)] = 0. (2.3) 

The LvN equation has been used to find the density operator for equilibrium systems that are 
stationary. Now the nonequilibrium system with their time-dependent coupling parameters 
follows the same equation, so the density operator (|2.1| ) that is directly defined in terms 
of the Hamiltonian itself does not necessarily satisfy the equation. In Sec. Ill we shall see 
how much the density operator satisfying Eq. ( |2.3| ) deviates from the instantaneous density 
operator ( |2.1| ) for time-dependent harmonic oscillators. 

In the context of quantum mechanics a powerful canonical method was developed by 
Lewis and Riesenfeld [p6| . They observed that any operator 0(t) satisfying the quantum 
LvN equation 

ihld(t) + [6(t),H(t)] = o, (2.4) 

can also be used to find the exact quantum states of Eq. (|2.2|) . In fact, any exact quantum 
state is given by 

\y(t)) = Y,Cne^ {t) \\ n ,t). (2.5) 

n 

where 

6(t)\\ n ,t) = X n \X n ,t), 

ln (t) = J dt(X n , t\ - H(t))\X n , t). (2.6) 

Another useful property of LvN approach is the linearity of the LvN equation. In fact, a 
product 6 1 {t)6 2 {t) satisfies Eq. whenever 6i(t) and 6 2 {t) satisfy Eq. Q. Therefore 
it holds that any analytic functional F[0(t)] satisfies Eq. ( |2.4j ) provided that 0(t) satisfies 
the same equation. In particular, we can still use 0{t) to define the density operator for the 
time-dependent system 

Po(t) = e —--, Z (t)=Tr[e- f}0{ % (2.7) 
Here /3 is a free parameter and will be identified with the inverse temperature only for the 



equilibrium system, in which the Hamiltonian itself satisfies Eq. (|2.4j ) and is used to define 
the standard density operator (|2.1|). 
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III. TIME-DEPENDENT OSCILLATOR 



As a simple nonequilibrium system we begin with an ensemble of time-dependent har- 
monic oscillators. This system carries more meaning than being merely a toy model because 
most of systems with some exception such as the massless $ 4 -theory, either in equilibrium 
or in nonequilibrium, can be approximated by a quadratic Hamiltonian around zero-force 
points, stable or unstable. In the technical aspect the time-dependent oscillator can be ex- 
actly solved in terms of its classical solution. We now focus on a general oscillator with a 
time-dependent mass and frequency squared 

^-tkf+^M' (3 ' 1) 

where uj 2 (t) is allowed to change the sign during a phase transition. The LvN approach will 
be employed below to find the exact Fock space and to construct the coherent, thermal and 
coherent-thermal states. 



A. Fock Space 

The key idea of the LvN approach to the Hamiltonian ( |3.1| ) is to require a pair of operators 

13,111 



a(t) = i(u*(t)p - m(t)u*(t)q\ 

a) (t) = -i(u(t)p - m(t)u(t)q) , (3.2) 



to satisfy the LvN equation ( |2.4| ). This results in the classical equation of motion for the 
complex auxiliary variable u 

u(t) + —rliiit) + uj 2 (t)u(t) = 0. (3.3) 
m{t) 

Note that these operators depend explicitly on time through u(t) and are hermitian conjugate 
to each other. Further, these operators can be made the annihilation and creation operators 
with the standard commutation relation for all the times 

[a(t),tf(t)] = l. (3.4) 

The above commutation relation is guaranteed by the wronskian condition 

hm{t)(u*{t)u{t) - ii(t)u*{tf) = i. (3.5) 

A comment is in order. Nothing prevents one from using these operators for an inverted 
time-dependent oscillator as far as Eq. ( |3.5| ) is satisfied. The inverted oscillator will be 
treated in detail in Sec. IV in the context of second order phase transitions. 

^From the argument in Sec. II, using a(t) and aJ(t) one may construct two particular 
operators that also satisfy Eq. (|2.4j ): the number and the density operator. By defining the 
number operator in the usual way 
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N(t) = tf(t)a(t), 

one finds the Fock space consisting of the time-dependent number states 

N(t)\n,t) = n\n,t). 



(3.6) 



(3.7) 



The vacuum state is the one that is annihilated by a(t) and the nth number state is obtained 
by applying cft(t) n-times to the vacuum state: 



a(t)\0,t) = 0, 

< 8, < ( »"|o, ( ). 



\n,t) 



In the coordinate representation the number state is given by (see Appendix A) 



(3.8) 



1 



1/4 



( u(t) ■ 

2irh 2 u*(t)u(t) ) ^n~\ \u*{t) 



2h 2 u*(t)u(t), 



exp 



i mil* (t) 2 



(3.9) 



where H n is the Hermite polynomial. 

Equation ( |3.2|) can be inverted to yield the position and momentum operators 

q = h(u(t)a(t) +w*(t)a t 
p = hm(t) (u(t)a(t) + u*(t)a f ) . (3.10) 

Hence the expectation value of the position and momentum with respect to each number 
state vanishes 



(n,t\q\n,t) = (n,t\p\n,t) = 0. 



(3.11) 



The only nonvanishing expectation values come from even powers of the position or momen- 
tum. The quadratic power of the position and momentum has the expectation values 



(n, t\q 2 \n, t) = h 2 u*(t)u(t)(2n + 1), 
(n,t\p 2 \n,t) = h 2 m 2 (t)u*{t)u{t){2n + l), 
(n,t\(qp + pq)\n,t) = h 2 m(t)(u*(t)u(t) + u*{t)u{t)){2n + 1). 

The Hamiltonian thus has the expectation value 

h 2 r i 
H n (t) = (n,t\H(t)\n,t) = —m(t) ii* (t)ii(t) + lu 2 (t)u* (t)u(t) (2n + V 



(3.12) 



(3.13) 



The prominent advantages of using a^(t) and a(t) appear when one tries, by using the 
standard technique of quantum mechanics, to construct other quantum states such as the 
coherent, thermal and coherent-thermal states. 
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B. Coherent State 



The coherent state is a particularly useful quantum state in quantum field theory for 
phase transitions or nonequilibrium dynamics. It may be used to obtain the effective po- 
tential for the classical background as on order parameter with contributions from quantum 
fluctuations. By treating a(t) and a'(t) as the annihilation and creation operators, we follow 
the definition of the coherent state for the time-independent oscillator [37]. The coherent 
state is defined as an eigenstate of a(t): 



a(t)\a, t) = a\a } t). 



(3.14) 



where a is a complex constant. The coherent state can be treated algebraically by intro- 
ducing a displacement operator 



D(a) 



-aa* (i)+a* a(t) 



which is unitary 



D(a)b\a) = D\a)D{a) = I. 



(3.15) 



(3.16) 



Moreover, the displacement operator translates by constants the annihilation and creation 
operators through the unitary transformation 



b(a)a(t)D ] (a) = a(t) + a, D(a)a ] \t)D ] (a) = a(t) + a*, 
and the inverse unitary transformation 

b\a)a{t)D(a) = a(t) - a, L> f (a) o f (t)D(a) = a(t) - a*. 



(3.17) 



(3.18) 



Hence one sees from Eq. ( 3.17 ) that the coherent state results from the unitary transforma- 
tion of the vacuum state 



\a,t) = D*(a)\0,t) 



a a x — A tx 

2 



(3.19) 



It has been known since Schrodinger that the coherent state gives rise to a classical field 
for the time-independent oscillator p8| . In our time- dependent case, from the expectation 
values of the position and momentum with respect to the coherent state 



q c (t) = (a, t\q\a, t) = h(a>u(t) + a*u*(t)^j , 
p c (t) = (a,t\p\a,t) = hm(t)(a>u(t) +a*u*(t) 



(3.20) 



one sees that q c indeed satisfies the classical equation of motion ( p.3| ), because u{i) has 
already satisfied Eq. (|3.3| ) and a is a constant. Besides, q c is real and p c = m(t)q c , so one 
may identify q c and p c with the classical position and momentum. 

In the above the coherent state has been constructed from the exact Fock space. There 
is another method to find the coherent state based on the minimization of action |Tj|. In 
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contrast with the LvN approach in which the operators ( |3.2|) are required to satisfy the 
LvN equation, one regards u as a free parameter, works on the ^-parameter Fock space and 
minimizes the Hamiltonian expectation value with respect to the coherent state. To show 
the method in detail, take the Hamiltonian expectation value with respect to the coherent 
state 



H c (t) = (a,t\H(t)\a,t) = H c (t)+H v (t), 
which consists of the classical part 

H c (t) 



1 9 m(t) 9, . 9 



2m{ty c ' 2 

and the vacuum fluctuation part given in Eq. (gig) with n = 

h 2 r 

H v (t) = —m(t)[u\t)u(t)+u 2 (t)u*{t)u(t) 
By writing the complex u in a polar form 

u (t) = ^l e -m 



(3.21) 



(3.22) 



(3.23) 



(3.24) 



in terms of which Eq. ( |3.5| ) becomes 9 = l/(2m^ 2 ), and by introducing = m(t)£, one 
obtains the effective Hamiltonian 



Hc(t) 



2m(t) 



Vl + ^\t)q 2 c + h 



1 



2m(t) 



P 2 + ^ 2 {t)e + 



i 



(3.25) 



The last term in the square bracket of Eq. ( |3.25| ) has the same form as the angular momen- 
tum for a particle having rotational symmetry in two dimensions, but its origin is rooted 
on the condition (|3.5| ) from quantization. Thus the effective Hamiltonian from the coher- 
ent state is equivalent to a two-dimensional Hamiltonian, which consists of the classical 
and quantum fluctuation parts. The variables q c and £ are independent and the Hamilton 
equations are for q c 



dq c 
dt 

dp c 
dt 



OHr 



dp c m(t) 
dH c 



Pc 



dq c 



m{t)u (t)q c 



(3.26) 



and for £ 



dt 
dt 



d{H c /h) _ 1 
dp^ m{t) 
d(H c /h) 



m(t)u/(t)£ + 



1 



(3.27) 



It is then easy to show that the Hamilton equations (|3.27|) equal to the second order equation 
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0. 



(3.28) 



and that Eq. (|3.28|) is nothing but the equation (|3~3|) when u has the form ( |3.24|) and satisfies 
the condition ( |3.5| ). Hence the minimization of the effective action gives the identical result 
as the LvN approach. Still another method is the mean field approach, in which the position 
and momentum are divided into a classical background and a fluctuation part 



Q = Qc + qf, P = Pc + Pf- 



(3.29) 



Then the total Hamiltonian is composed of three parts: the classical background and fluc- 
tuation parts 



H(t)=H c (t) + H f (t)+H int (t), 



where 



HAt) 



2m{t) 



p) + ^\t)q) 



(3.30) 



(3.31) 



is the fluctuation Hamiltonian, and 



Hi n t{t) 



m(t) 



p c p f + m(t)u 2 (t)q c q f 



(3.32) 



is the interaction Hamiltonian between the classical background and the fluctuation. And 
then quantize the fluctuation Hamiltonian ( 3.31 ) according to the method in the previous 



subsection but keep the classical one unquantized. Since the last two terms proportional to 
qj and pf have the zero expectation value, the expectation value of the total Hamiltonian 
( |3.30| ) with respect to the vacuum state of the fluctuation Hamiltonian (|3.31|) yields exactly 
the effective Hamiltonian ( |3.25| ). Therefore, it has been shown that the expectation value of 
the original Hamiltonian with respect to the coherent state is equivalent to the sum of the 
classical part fl3.22p and the vacuum expectation value of the fluctuation part (|3.31|). 



C. Thermal State and Density Matrix 

The ensemble of time- dependent oscillators exhibits intrinsically nonequilibrium behav- 
iors, so it does lose a rigorous physical meaning to attribute any thermal property to the 
density operator fl2.1|) . However, the harmonic oscillator problem is exactly solvable, so even 
in the time-dependent oscillator case one may look for a density operator which is quadratic 
in the position and momentum, and then fix their variable coefficients to satisfy the LvN 
equation |)39|| . On the other hand, in the LvN approach we can still use the operators 



that have already satisfied the LvN equation and define the density operator from them. But 
there still remains a free parameter to incorporate the initial thermal equilibrium. We study 
the physical meaning of the density operator and see how the initial thermal equilibrium 
evolves quantum mechanically . 

By noting that N(t) satisfies Eq. ( |2.4|) , we define the density operator by 
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Pt 



Z N 



(3.33) 



where (3 and ujq are free parameters and is the partition function given by 



EMI 



-phw (N{t)+\) 



n,t) 



1 



n=0 



2sinh(^; 



(3.34) 



It has the same form as the standard density operator, the time-independent annihilation 
and creation operators now being replaced by the time-dependent ones fl3.2p. So Eq. fl3.33| ) 
includes the time-independent special case by choosing (3 = l/(k B T) and u the 

oscillator frequency. In the coordinate representation the density matrix is given by (see 
Appendix B) 



^ n (q',t)%(q,t)e 



tanh(^^ 1/2 
2nh 2 u*u 

x exp 



exp 



4hdt Huu)iq ~ q > 



8h 2 u*u 



(q' + q) 2 tanh(^) + (q> - qf coth(^; 



(3.35) 



Now the density matrix ( |3.35| ) can be compared with that for the time-independent oscillator 
and the density operator (|2.1|) for the instantaneous Hamiltonian. For that purpose we 
restrict our attention to the particular case, in which the mass is constant, m(t) = mo, 
and uJ 2 (t) is positive (see Sec. IV for the sign changing case of uJ 2 (t)) and slowly changing 
\ui(t) / uj(t) \ <C 1. In that case we may look for the solution to Eq. (|3.3|) of the form 



u{t) 



1 



2hm Q( y t) 



-i J n(t) 



(3.36) 



where 



n 2 (t) = u 2 {t) + 



3Q 2 (t) lQ(t) 



(3.37) 



The adiabatic (WKB) solution is obtained by approximating Q(t) « u(t). Then the density 
matrix ( |3.35| ) reduces to the adiabatic one 



pA(q',g,t) 



m uj(t) tanh(^a^ 1/2 



x exp 



rriQU}(t) 



exp 



4 h u{t) 



(q 



(q + q) 2 tanh( 



+ (q ~ q) 2 coth( 



/3hu , 



(3.38) 



On the other hand, the density operator ( |2.1| ) for the instantaneous Hamiltonian has the 
matrix representation 
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m to(t) tanh(^p^ 



x exp 



moUj(t) 
4/T~ 



1/2 



'(,' + g) 2 tanh(^l) + (,' - g) 2 coth(^; 



(3.39) 



In the case of the time-independent oscillator with u(t) = ujq, the density matrices ( |3.35| ) 
and ( |3.39| ) reduce further to the standard one |40| . The instantaneous density matrix is 
compared with the adiabatic one by taking the ratio 



PH(?',g,t) 

PA(q',q,t) 



tanh( 



2 - 



tanh( 



1/2 



exp 



x exp 



Ah 



(q + qf 



i m Co{t) 
sinh(f {u(t)-u 



(q 2 



cosh(- 



cosh( i 



- (q ~ qf 



sinh(f (u(t)-u ] 



sinh(- 



sinh( / 



(3.40) 



The second factor in Eq. ( |3.40| ) gives rise to a small phase factor because u(t) is slowly 
varying. As far as u(t) remains close to u>o, the instantaneous density matrix (|3.39| ) is close 
to the adiabatic one (|3.38|) . Otherwise, the exact nonequilibrium evolution (|3.38|) is far away 
from the quasi-equilibrium one described by ( p.39|) . This implies a significant deviation of 
the nonequilibrium evolution from the equilibrium one as u(t) differs from uj by a large 
amount. 

Being mostly interested in phase transitions, we assume the system to start from an 
initial thermal equilibrium at early times. This requires that the oscillator have a constant 
frequency ujq at early times and j3 be fixed to the inverse temperature. As in the case of co- 
herent state, the evolution of the initial thermal state can be found the effective Hamiltonian 
from this thermal state. From the expectation values (see the appendix of Ref. [p8R ) 



(P 2 h 



Tr 
Tr 



p T {t)q 2 ] =ftV(iMt)coth(^p), 
p T (t)p 2 } = h 2 m 2 (t)u*(t)u(t)coth(f^) 



(3.41) 



one obtains the effective Hamiltonian from the thermal state 



H T (t) = Tr p\(t)H(t) 



h 2 M ..ufP^o* 
—m[t) coth( — - — 



u*(t)u(t) + co 2 (t)u*(t)u(t) 



(3.42) 



Once again by using the complex parameter u in the polar form (|3.24|) and by introducing 
= m(t)(, one may rewrite the effective Hamiltonian as 



H T (t) = frcoth( 



2m{t) 



pI + ^u 2 (t)( 2 + 



8m(t)f 2 



coth(^^)F v (t). (3.43) 



Since hcoth.(j3hu)o/2) is constant, H-T(t) has the same Hamilton equations ( |3.28| ) as ifv(£)- 
This is identical to the classical equation of motion ( |3.3|) together with the boundary condi- 
tion (JjHfy. 
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D. Coherent-Thermal State 



A more general density operator that is at most quadratic in the position and momentum 
was introduced in Ref. p8fl . It has the form 



Pc.tO) 



1 



exp 



-(3{fkv tf(t)a(t) + 5a ] {t) + 5*a(t) + e } 



(3.44) 



where Zc.t is a partition function. In fact, the density operator of Eq. ( |3.44| ) can be 
transformed into that of Eq. (|3.33|) by the unitary transformation 



D\a)p c (t)D(a) = p(t), 



(3.45) 



where £>{a) with a = 5/{hu ) and e = (Huj /2) + (\5\ 2 /hu) ). From the expectation values 



2 )c.t = Tr 



'C.T 



Tr 



PCT(t)q 2 } = ql + h 2 u*(t)u(t) coth(^p), 
PCT(t)p 2 } =Pl + h 2 m 2 (t)u*(t)u(t) coth(^), 



(3.46) 



follows the effective Hamiltonian 



H c . T (t) = H c (t) + H T (t), 



(3.47) 



where H c is the classical Hamiltonian. The effective Hamiltonian ( |3.47|) has almost the 
same form as Eq. ( |3.25| ) except for the overall factor coth(/?^o/2), hence describes the 
same Hamilton equations (|3.27| ). 

In summary, what we have shown in this section is that the LvN approach provides 
us with the exact Fock space and various quantum states for time-dependent oscillators. 
The only auxiliary field necessary for that purpose is a complex solution u to the classical 
equation of motion (|3.3j ) that satisfies the wronskian condition ( |3.5| ) from quantization. It 



has further been shown that the LvN approach is equivalent to the minimization principle 
for the effective action and to the mean field method. However, the advantage of the LvN 
approach lies in the manifest correspondence with quantum mechanics and quantum many- 
particle system and in the readiness to apply their standard techniques. For instance, the 
density matrix (|3.35|) has been found, which has many features similar with the standard 



one. 



IV. INVERTED HARMONIC OSCILLATOR 

As a toy model for the second order phase transition, let us consider the time-dependent 
harmonic oscillator 

H = \f + ^ 2 (t)g 2 , (4.1) 

where w 2 (t) has the asymptotic value w 2 (> 0) far before and —uj 2 {< 0) far after the quench. 
The conspicuous point of the model is the sign change of uj 2 (i). At earlier times before the 
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quench the oscillator executes a stable motion about q = 0, the global minimum, but after 
the quench the potential is inverted and q = becomes an unstable configuration. 

At earlier times far before the phase transition, the complex solution to Eq. ( |3.3| ) satis- 
fying Eq. (|3.5| ) is given by 



Ui{t) 



-iuiit 



(4.2) 



According to Eq. ( |3.2|) , the Fock space is now constructed from the annihilation and creation 
operators 



a(t) 



D iuJit 



\JTiuJi 

^—iuJit 



(ip + Uiq 



e a , 



(-ip + uo iq ) = e-^al 



(4.3) 



Note that a(t) and a)(t) differ from the standard ones do and a\ only by phase factors. 
Though the Hamiltonian has the standard representation 



Hi = hu)i(a ] (t)d(t) + = huji(ala + ~ 



(4.4) 



the phase factors are necessary for a{t) and <${t) to satisfy the LvN equation. Hence the 
vacuum expectation value is given by the well-known result 



H 



v 



-Tiojj 



and the coherent state ( |3.19| ) yields 



2 i 2 2 , 1 * 



Now the density operator (|3.33|) reduces to the standard one 



Pi,T = 



-/3?kJi(aJa + ±) 



(4.5) 



(4.6) 



(4.7) 



after identifying ujq = Ui and (3 = l/iksT), and leads to the Hamiltonian expectation value 

(4.8) 

On the other hand, at later times far after the quench, the solution to Eq. fl3.3p is given 



rr I* ^{P huJ i 

^i,T = -dWjCOth^ — - — 



by 



«/(*) 



2h 
l 



Ci(u>i, ujf) cosh(cj/t) — iC^LOi, ujf) sinh(cj/t) 
(C^Uf) -iC^u^e^ + (diuuUf) +iC 2 {u i ,u f ))e~"f t ], (4.9) 
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where Cj,j = 1,2 depend on the intermediate process toward the final state. Remarkably, 
the vacuum and thermal expectation values vanish: 



H 



H 



0. 



(4.10) 



since the kinetic and potential energies contribute equally. Whereas the expectation value 
with respect to the coherent state ( |3.19| ) and the general density operator (|3.44| ) takes the 
form 



H 



f,c 



H 



f,C.T 



-0J 



• 2 n 2 



(4.11) 



This implies physically that as the system undergoes the phase transition out of equilibrium, 
quantum effects vanish and it becomes classical. Large uncertainty has been suggested as 
one of the criteria on the classicality 



n 



(Aq)(Ap) = — ujfC l (oji, io f) cosh (ujft) + ujfC l (cUj, uf) sinh (ujft) 



(4.12) 



Though the oscillator starts with the minimum uncertainty given by Eq. ( [4 . 2| ) , its uncer- 
tainty increases exponentially and it becomes eventually classical after the completion of 
quench. When Eq. ( P~9| ) is substituted into Eq. ( |3.35| ), the density matrix depends on the 
intermediate process and spreads as ^Ju*j(t)uf(t). Therefore the quintessence of second order 
phase transitions lies in the whole process how systems evolve out of equilibrium from their 
initial equilibrium. 

To show explicitly how the nonequilibrium dynamics depends on the intermediate pro- 
cesses, we consider two exactly solvable models. The first model, which is the zero mode of 
the free scalar model in Sec. VI, is an oscillator describing a finite smooth quench with the 

mass 



m 2 (t) 



t 



Iq tanh(— ) . 



(4.13) 



The mass has m 2 = + m\ at earlier times (t = — oo) and has — = — TOq + m\ < 
at later times (t = oo). r measures the quench rate, i.e. the rate of change of mass. The 
instantaneous quench corresponds to the (r = 0)-limit. The solution to Eq. ( |3.3| ) is found 



u{t) 



-m.it 



m f),~2( im i 



m f )]l 



D 2t/r 



^JITim^ ± ^ 2 

At earlier times the solution (f4.14|) has the correct asymptotic form 

p-im.it 

Ui(t) 



On the other hand, at later times the asymptotic form of Eq. ([4.14|) becomes |42| 

1 (-l)r(l - imiT)T(m f r) 



(4.14) 



(4.15) 



yj2hmi \ {imi — mf)T 2 



e m f t 



-^{irrii - rrif) 
1 (-l)r(l - imiT)r(m f r) 
y / 2hm i ^(irrii + m,f)T 2 (— \{imi + m/) 



-met 



(4.16) 
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Two points are observed: the initial solution branches into an unstable growing and a 
decaying mode as expected and the coefficients C\, C 2 of Eq. ( |4.9|) depend on the the mass 
parameters m^mj and the quench rate r. In other words, the final asymptotic state of 
nonequilibrium evolution depends on the intermediate process. 

The next model describes various quench processes and exhibits how the stable mode 
far before the quench branches into the unstable growing and decaying modes during the 
quench. Without loss of generality the nonequilibrium phase transition is assumed to take 
place through the time-dependent frequency (mass) squared 



u\t) 



"I 



(2Zi+l)/(2Z 2 +l) 



UJ 



2| t -t 
to—ti 



(2Zi+l)/(2Z 2 +l) 



—UJ 



2| t-tp 
tg—ti 



ti > t. 



to > t > ti, 



tf>t> t , 



(4.17) 



, ,2 — , ,2 1 tf-tp 



(2Zi+l)/(2/ 2 +l) 



t > t 



where l\ and l 2 are non-negative integers. Here to — ti adjusts the rate of and tf — to the 
duration of the quench. The particular case of l\ = l 2 = is used as the finite linear quench 
model [13,16|. Before the time to, the system maintains the symmetry about q = 0, the 
minimum of the potential. But as time goes on after to, q = remains no longer the true 
minimum of the system and the symmetry is broken. The particular form of the power-law 
in Eq. (|4.17|) is chosen to allow an analytical continuation of oo 2 (t) for changing the sign and 
to make its derivatives also continuous. Before t, > t, the solution is given by Eq. ( fOj) . 
During to > t > ti, the solution is given by a linear superposition of Hankel functions 



u(t) = D x z v Hf\z) + D 2 z v Hj / 1 \z), 



where 



tp-t 
to — ti 



l/2v 



2/ 2 + l 



2h + Al 2 + 3 



(4.18) 



(4.19) 



Here and are positive and negative frequency solutions, respectively. The constants 
Di and D 2 are determined by continuity of u(t) and it(t) across tf 



v / 2hu i Az L i 



d 



-iz, 



2(t - U) 



^JZhhJi Az\ 



A 

dzi 



Hi 2 \z i ) + {uj l 



2(t - U) 



(4.20) 



where 



Zi = 2uui(t - U). 



(4.21) 



Beyond the quench time to, it is necessary to do carefully the analytic continuation and to 
take a suitable Riemann sheet so that oj 2 (t) and its derivatives are to be continuous from 



> to uj 2 < across to- The analytic continuation of the solution ( [4.1 8|) yields 
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u (t) = 



(4.22) 



where 



t-t 
to — U 



l/2u 



(4.23) 



At later times (z ^> 1) during the quench, the solution (|4.22 ) has the asymptotic form [43 



Die 



(4.24) 



Thus the stable mode of the oscillating solution ( |4.2j ) branches into the growing mode ( |4.24| ) 
which dominates during various quench processes and into the decaying mode which con- 
tributes negligibly to the correlation functions. The asymptotic solution (|4.24|) also depends 
on the intermediate processes through t ,ti,tf and h,h- 



V. TIME-DEPENDENT ANHARMONIC OSCILLATOR 

In this section we extend the formalism developed in Sec. Ill to the time-dependent 
anharmonic oscillator with the Hamiltonian 

2 

H(t) = ^r + m(t)V(q), (5.1) 

where 

n?) = ^f§? 2 " (5-2) 

Though the potential of a power law is assumed, the formalism can readily be generalized to 
any polynomial and analytic potential. In the time- independent case (\2 n = constant), the 
variational perturbation method has been introduced as one of the powerful methods to find 
the Hilbert space ||44j| . The vacuum state in this approach is the Gaussian wave functional 
that minimizes the effective action. The excited states are then obtained from the vacuum 
state just as number states of a harmonic oscillator are obtained from the ground state. 
Though these states can be calculated explicitly in terms a complex solution to the classical 
equation for motion, they are in fact equivalent to those from the mean-field method. 



However, there have been some attempts to go beyond the Gaussian state. In Ref. |45 
a scheme was proposed to find the operators, which generalize the annihilation and creation 
operators and satisfy the LvN equation (|2.4j) , to all the orders of coupling constant in the 
time-independent case and to the first order in the time-dependent case. In particular, 
the generalized annihilation and creation operators for the time-independent oscillator with 
a quartic potential satisfy a q-deformed algebra rather than the standard commutation 
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relation ||46|| , from which follows an algebraic construction of excited states and energy 



spectra beyond the variational Gaussian approximation. It would be interesting to find 
such an algebraic structure for interacting quantum fields, which may shed some light on 
the nonperturbative method beyond the mean-field method. Also it would be interesting 
to compare this scheme with other nonperturbative methods in the time-independent case 
such as the perturbative expansion method around the Gaussian effective action [|47]] and 
the time-dependent variational method [JS|. But we shall not pursue further this issue in 
this paper. 



A. Fock Space 

In the case of time-dependent anharmonic oscillators, the LvN approach searches for the 
annihilation and creation operators that are still linear in the position and momentum and 
satisfy the LvN equation 



A(t) = i(y*(t)p - m(t)v*(t)qj, 

A ] {t) = -i(v{t)p - m(t)v(t)qj. (5.3) 
One then requires them to satisfy the LvN equation fl2.4|), leading to the equation 



i>m+ m im+v{t) s m =0t (5 . 4) 

and further differentiates functionally with respect to q and takes the vacuum expectation 
value of the resultant equation 

v(t) + ^v(t) + (0, t \ 6 -^-\0,t)v(t) = 0. (5.5) 

Here the vacuum state is annihilated by A(t) 

A(t)\0,t)=0. (5.6) 

One makes A(t) and A^(t) the annihilation and creation operators, respectively, by imposing 
the standard commutation relation for all times 

[A(t),A\t)] = l. (5.7) 

Equation (|5.7|) is equivalent to the wronskian condition 

hm{t)(v*{t)v{t) -v{t)v*{t)) =i. (5.8) 

The Fock space consists of the number state obtained by applying A\t) n-times to the 
vacuum state 

\n,t) = V '> |0,t). (5.9) 
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These number states are excited states and have the coordinate representation ( |3.9| ) now 
with u(t) replaced by v(t). 

^From the position and momentum operators expressed in terms of the annihilation and 
creation operators 



(5.10) 



q = %(y(t)A(t) + v*(t)A^(t)), 
p = Tim{t)(y(t)A(t) + i)*{t)A\€j), 
follow the vacuum expectation values 

(2n)! 



(0,t\q 2n \Q,t) 



2 n n\ 



H 2 v*(t)v(t) 



(0,t\p 2 \0,t) = h 2 m 2 (t)v*(t)v(t). 
Then the classical equation of motion ( |5.5|) becomes 



V{t) + ZT7I\ v (' t ) + 



ra-1 



h 2 v*{t)v{t) v{t) = 0. 



m{t) 2 n ~ 1 (n - 1)! 

There is another method to derive Eq. (|5 . 1 2|) . By using the Wick-ordering 

(2n)\Ti 2n 



~2n 



q 



E 



v*(t)v(t)] k : [v(t)A(t) + v*(t)A\t) 



2{n-k) 



k ^2 k k\{2n - 2k)\ 

one obtains the Hamiltonian truncated at the quadratic order of A(t) and A(t) 

jj _ J r.-..,i4\ t;,i-i-\ -i;, ii\;,i4\ v: n- \ a n- \ • ■ ■ ■>■•-. IT/+^^2" 



/ ( ; = -/i 2 m(t)[(t;(t)i(t)) 2 + 2 ? )*(t)t;(t)it(t)i(t) + (^(t)it( t ))^ 



+ m(t) 



h 2n X 2n (t) r ..„ 

2"(n - 1)! 



(5.11) 



(5.12) 



(5.13) 



(t)^(t)]"- 1 [(^(t)i(t)) 2 + 2v*{t)v{t)A\t)A{t) + (t;*(t)it(t))2] j ( 5 . 14 ) 



where purely c-number terms are neglected. One then requires A\t) and A(t) to satisfy the 
LvN equation ( |2.4j ) for the truncated Hamiltonian Now the LvN equations for A'(t) 

and A(t) lead exactly to the equation of motion (|5.12|) . 



B. Effective Hamiltonians 



Still another method to derive the equations of motion for v and q c is the minimization 
principle for the effective action. For that purpose we consider a complex v (t)-parameter 



family of the Fock spaces constructed by the annihilation and creation operators ( |5.3|) and 
take the Hamiltonian expectation value with respect to various quantum states such as the 
vacuum, coherent, thermal and coherent-thermal states. We do not require A{t) and A\t) 
to satisfy the LvN equation a priori, but minimize the action to determine the equation of 
motion for v(t). 

The first effective Hamiltonian is the vacuum expectation value 



H Y (t) = (0,t\H(t)\0,t) = ^m(t)v*v + m(t)^ 



h 2 v*(t)v(t) 



(5.15) 
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The next state under consideration is the coherent state, which is obtained by applying a 
displacement operator to the vacuum state, 

\a,t) = &(a)\0,t) = e ait(t) - a * iw |0,t). (5.16) 

Then the coherent state expectation value leads to the effective Hamiltonian 

H c (t) = (a,t\H(t)\a,t), (5.17) 

which, with the aid of Eq. ( |5.13| ), is decomposed into 

H c (t)=H c (t) + H q (t). (5.18) 

Here H c (t) is the classical part 



(5.19) 



2m{t) v ' (2n) 

and H q (t) denotes all the quantum contributions including the vacuum expectation value 

fl 2 , A 2n (^) r. o , , i k 



H,(t) = - m (t)i'i + „(«) g 2U!(2 ;^ M)[ [ftV(«Mt) 
where 

g c = (a,t|g|a,t), p c = (a, t|]3|a, t). 
The final state is the thermal state defined by the density operator 



q, 



2(n-k) 



PT 



1 c -0u, n (AHt)A(t)+±) 



J N 



(5.20) 



(5.21) 



(5.22) 



with Zn being the partition function. The density operator (|5.22 ) leads to the effective 
Hamiltonian 



h 



H T (t) = Tr p T H(t) = -m(t)v*v + m {t)^-{q 2 ) 



2 n n\ 



where 



(g 2 ) T = ^V^COth(^). 



(5.23) 



(5.24) 



Likewise, the density operator of the form (3.44) with a)(t) and a(t) replaced by A^(t) and 
A(t) leads to the effective Hamiltonian 



ti 2 



H c . T (t) = H c (t) + —m(t) coth( o 



IV V 



\ 2 k k\{2n-2k)\ 



[q 2 )Wc in - k) - 



(5.25) 



Note that Eqs. (|5.15| ) and ( |5.18| ) are also obtained from Eqs. ( [5.23] ) and (|5.25| ), respectively, 
by replacing h 2 v *v coth(/?a;o/2) with h 2 v*v or taking the zero-temperature limit (j3 — > oo). 
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We now study the dynamics of the effective Hamiltonians. We mainly focus on the 
effective Hamiltonian ( 5.25 ) since Eq. Q5.18 ) is the limiting case of Eq. ( 5.25Q when f3 — > oo, 
i.e., T — > 0, and Eq. ( p,15| ) is the limiting case of Eq. ( |5.1S| ) when q c = p c = 0. The 
independent variables of the Hamiltonian (|5.25| ) are (q c ,p c ), (v,p v = mi)*) and (v* ,p* = mv). 
So we obtain the equation of motion for q c 



m 



m {2n — 1)1 



2n-l 



4- V ^n{t) ( 2 \k 2n-2fc-l _ () 

+ ^W.(2n-2k-l)\ {9 Mc 



(5.26) 



The equation of motion for v* is given by 



III 



{2 k - 1 (k-l)\{2n-2k) 



2n-2k I ^2\fe-l 



<r>?> = o, 



(5.27) 



and the complex conjugate of Eq. (|5.27| ) is for v. The equations of motion from the effective 
Hamiltonian (|5.18 ) is the limiting case of Eqs. ( |5.26| ) and (|5.27 ) when (<? 2 )(t) = fi 2 v*v: 



m 



q c H — Qc + 



A 2 „(*) 



,2n-l 



m (2n — 1)! 



+ E 



2 *„,\fc„2n-2fc-l 



\ 2 k k\(2n -2k- 1)! 



(h z v*v) k q i 



m . x - 

f H f + > —, — r 



A 2 „(t) 



gf - 2fe (^^) fe -^ = 0. 



r;2 fc " 1 (A;-l)!(2n-2A;)! 



(5.28) 



And the limiting case q c = of Eq. (|5.28|) is the equation for the effective Hamiltonian 

(EH) 



A2n(*) 



m 

m 2"- 1 n-l)! 



(H 2 v*v) n - L v = 0. 



(5.29) 



Note that Eq. Q5.29| ) is identical to Eq. ( |5.12j ) from the LvN approach. 
Or, by writing v in the polar form 

v(t) = ^e-»®, 



(5.30) 



and by introducing the momentum = m(t)(, the effective Hamiltonian (|5.25|) is rewritten 

as 



Pc 



2m(t) 



(2n)\ 

n 



p\ 



2m(t) 8m(t)C 2 



A2n(^) 



\2 k k\{2n - 2k)\ 



h( 2 coth( 



c 2(n-k) 
He 



(5.31) 



Now the equation of motion for the classical field q c is given by 

) 



m 



(2n - 1)! 
and that for £ by 



0, (5.32) 
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< + £ - W + £ ^- 1 ( fc -l)!(l-2 fe )r ^- 2t ^ 2 «W£)]"C = 0. (5.33) 
The phase of v (t) is obtained by the integration 

w = f*zmry (5 ' 34) 



C. Coherent State vs. Hartree-Fock Method 



In this subsection we show that the nonequilibrium dynamics obtained from the effective 
Hamiltonian ( |5.25[ ) in the LvN approach recovers exactly the equations of motion from the 
mean field and Hartree-Fock methods. 

First, the effective Hamiltonian from the coherent state can be obtained in another way. 
By dividing q and p into a classical background and a quantum fluctuation 



Q = Qc + Qf, 



(5.35) 



we obtain the expectation value for the Hamiltonian ( |5.1| ) with respect to the thermal state 
(PI) 



2m{t) 2m(t) 
where we have used 

n 

(? n h = E 

and 



(V 2 f )T 



A 2 n(t) 



2 k k\(2n-2ky w/ 



2\fe 2(n-k) 



mm 



(2n)\ 



k=0 



2 k k\(2n - 2k)\ 



2(n-k) 



% v*v coth( 



(3hu> , 



(qf +i) T = = (pf +1 ) T . 
Noting that k = term recovers the classical potential and 

/~2\ 2+2 •*• ,u/^ W 0\ 

{Pf/T = m n v wcoth( ), 



(5.36) 



(5.37) 



(5.38) 



(5.39) 



we can show that Eq. (|5.36| ) coincides with Eq. ( 5.25 ). Therefore the coherent state leads 
exactly to the result from the mean-field method. 

Second, the Hartree-Fock factorization theorem ffJIJl leads to the effective Hamiltonian 



#h.f(*) 



"W /H.F) 



2m(t) 



+ m(t)- 



(2n)\ 



(5.40) 



where 



(P 2 }r.f = pI+P}, 

(g 2 )H.F = q 2 + q 2 f 



(5.41) 
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and 



(2n)! 



2n-k~k 

Q c Qf 



^ Q 2 k k\(2n - 2jfe)l 



g 2(n-fc) 



(2n)! 



n-l 

+ 2 ¥k\(2n - 2k - 1)! 



2n-2fc-l /~2\fc ~ 
<?c {Qf/TQf- 



for n > 2. The thermal expectation value of the equation of motion for q c yields 



m . -A A 2n (i) / A 2\fc„2n-2fc-l 



g c + — q c + 2^ 



k=0 



and for g/ 



g/ + m g/ + ^ 2 fc " 1 (A;- l)!(2ra - 2fc)! W/ 



2 fc fc!(2n - 2fc)! 



A 2 n(t) 



2 ) ^l g 2(n- fe )^ = q 



fc=0 



(5.42) 



(5.43) 



(5.44) 



Therefore it has been shown that Eqs. ( |5.43| ) and Q5.44| ) are the same as Eqs. ( |5.26| ) and 
( |5.27| ) from the coherent state representation. 



VI. FREE SCALAR FIELD FOR PHASE TRANSITION 



As a simple field model for the second order phase transition, we consider a free complex 
scalar field, the mass of which changes the sign during the quench.[] The system is described 
by the Lagrangian density 



£(x,t) = $*(x,t)$(x,t) - V$*(x,t) • V$(x,t) -m 2 (t)$*(x,t)$(x,t). 



(6.1) 



Here the coupling parameter m 2 (t) is assumed to begin with an initial positive value before, 
to change the sign during, and to reach a final negative value after the quench. The $ and 
$* are treated as independent fields. The Hamiltonian is given by 



H(t) 



d 3 x 



where 



n*(x, t)II(x, t) + V$*(x, t) ■ V$(x, t) + m 2 (t)$*(x, t)$(x, t) 



n(x,t) = ^4 = $*(x,t), 



(6.2) 



a$(x,t) 
n*(x,t) = i^ = $(x,t) 



(6.3) 



x The complex scalar field model may be related with the time-dependent Landau-Ginsburg theory 
provided that the free energy be interpreted as the Hamiltonian in this paper 
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are conjugate momenta. 

The field and momentum are Fourier-decomposed as 



*(x,t) 



d 3 k 



(2tt) 3 

n(x,t) = $*(x,t) 



0k(t)e lk - x , 



^ 3 k • 



and their conjugates as 



$*(x,t) 



d 3 k 



(2tt) 3 

n*(x,t) = 6(x,t) = 



(2tt 



rf 3 k 

(2^ 



vr k (t)e" 



-ik-x 



(6.4) 



rf 3 k 
(2^ 



k (t)e 



ikx 



rf 3 k 

(27rp 



ikx 



(6.5) 



So space integrals of quadratic fields and momenta result in momentum integrals for the 
decoupled modes 



J d 3 x$*$ 

d 3 xirn 



d 3 k * 
(2^ 0k0k ' 
ci 3 k 



cf xV$* ■ V$ 



(2tt) 3 
rf 3 k 
(2T) 3 



k 2 0*0 k . 



(6.6) 



One then obtains the Hamiltonian as the sum of infinite number of time-dependent harmonic 
oscillators 



H{t) 



(X) 3 



7T k 7T k + (k 2 +m 2 (t))0 k k 



(6.7) 



Canonical quantization is prescribed by imposing the commutation relations at equal times 

[</> k '(t),7T k (t)] = iM kjk /, 
[0 k ,(t),7T k (t)] = ih8 k y, 



(6.8) 



and all the other commutators vanish. Following Sec. Ill, we find the two pairs of the 
annihilation and creation operators ( |3.2| ) for each k-mode, 



a k (t) = i(<Pk(t)K ~ 0k(*A 



and 



c(*) = -*(Vk(t)vTk - ^k(t)&), 



«kW = <^(^k-0kW0 k ), 



(6.9) 



(6.10) 
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where (p k and tp k satisfy the same classical equation of motion 

k (t) + (k 2 + m 2 (t))ip k (t) = 0. 
They further satisfy the standard commutation relations 

[ak',a k ] = <5 k)k /, [a k /,a k ] = 5 k)k '. 



(6.11) 



(6.12) 



The Fock space for each mode can be constructed according to Sec. III. We consider two 
symmetric states: the vacuum and thermal states. The vacuum is the one annihilated by 
all the a k and a k : 



a k (t)|0,t) =0, at(t)\0,t) = 0. 
By inverting Eqs. (|6.9| ) and (|6.10|) one expresses the fields as 

0k = ^-(^kflk + y?k*a k V 
k = h(<Pka*k + y^k*a k y 
from which follow the vacuum expectation values 

d 3 k r 



(6.13) 



(6.14) 



(n*n) 



V 



(2tt) 3 
(2^)3 



(6.15) 



The initial thermal state defined by the density operator for each mode 

1 



p(t) = UM t ) = U ^ ex P 



1 

x — exp 



(6.16) 



leads to the thermal expectation values 



($*$) T = Tr 

(n*n) T = Tr 



p(t)$*$ 

p(t)fl*fl 



d 3 k 

(2^)3 
(27p 



h'<p* k (t)<p k (t) coth( 



i,k n 



K<pl(t)<p k {t) coth( 



2 ' 

(3huj i)k , 



(6.17) 



One then finds the two-point correlation functions at equal times by taking the expectation 
value with respect to the vacuum state 

d 3 k 



Gv(y,x,t) = ($*(y,t)$(x,t)) v 
and with respect to the thermal state 

G T (y,x,t) = ($*(y,£)$(x,t)) T = 
where 



(27T) 



3 tk-(x-y) 



^V k (t)^ k (t)coth(^; 



(2tt 



3 ik-(x-y) 



w i( k = V k2 + ^ 2 (-oo)- 



(6.18) 



(6.19) 



(6.20) 
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A. Instantaneous Quench 



The instantaneous quench model is an analytically solvable one, in which the mass 
changes as 



m 2 (t) 



-mj, t > 0. 



(6.21) 



Before the quench (t < 0), the solution to Eq. (|6.11|) , which also satisfies the condition ([T| 
is given by 



-e~ l ^\ uj L]l = \/k 2 + m 2 



(6.22) 



i.k 



Then the two-point vacuum correlation function ( |6.18 ) becomes 

d 3 k h 



G f i ,v(y,x,t) = J 



3 ik-(x-y) 



( 27r ) 3 2yjk 2 + mf 
ft miiTi(mi|x-y|) 



4tt 2 



l x ~ y 



(6.23) 



where K\ is the modified Bessel function. Equation ( |6.23|) coincides with the result for a 
massive scalar field in Ref. [HUJ. Similarly, the two-point thermal correlation function is 
given by 



G iiT (y,x,t) = G iiV (y,x,t) 

tl ITLj 



+ 



2tt 2 



x - y| + m i n=i 



^^KJlx-y^ + ^n) 2 ))- 



(6.24) 



On the other hand, after the quench (t > 0), the classical equations of motion are 
classified into two types: the one from the long wavelength modes with k 2 < m 2 has the 
negative frequency squared and exhibits an exponential behavior, and the other from the 
short wavelength modes with k 2 > m 2 still has the positive frequency squared and shows 
an oscillatory behavior. Each mode moves under a constant frequency squared before the 
quench time but suddenly experiences a potential step in the case of short wavelengths and 
a potential barrier in the case of long wavelengths. There is the analogy between Eq. (|6.11|) 
and the scattering problem of quantum mechanics. The solution to Eq. ( |6.11| ) together with 



the initial asymptotic data (|6.22|) is the complex conjugate of the scattering wave function 
by either the potential step or barrier |56| . The solution to Eq. ( |6.11| ) after the quench 
should match at the quench time t = continuously with Eq. ( |6.22|) before the quench. It 



is rather straightforward to find such solutions for the short wavelength modes (k > m . 



ffsM 



k / \ / \ 

-i — — sm(u;j k t) + cos(u;/k^) 



w /,k 



k 2 — m 2 , 



(6.25) 



and for the long wavelength modes (k 2 < in 2 . 
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i.k 



—l- 



sinh(o; f^t) + cosh(u f^t) 



mj — k 2 . 



(6.26) 



A few comments are in order. The solution (|6.26|) represents an instability due to the 
phase transition and is obtained by continuing analytically the solution (|6.25|) . When Eq. 
( |6.25| ) is rewritten as 



V/,k(*) 



2fuu, 



.k 



2w 



2cu 



/■k 



(6.27) 



the first and second terms correspond to the positive and negative frequencies, respectively, 
hence the second term explains the particle creation by changing frequency [[3^] . 

After some manipulations of algebra, we obtain the two-point vacuum correlation func- 
tion after the quench 



G/,v(y, x, t) = G it v(y, x, t) 

h r m f 



47r 2 |x — y| Jo 

h 



dkk 



/,k 



sin (A; | x — y|) sinh 2 (uj^t) 



i ~ 2 

\ ^t.k^k 



. sin(A;|x - y|)sin (co ftk t). 



47r 2 |x — y 

Similarly, the two-point thermal correlation function is given by 



(6.28) 



G /,T(y,x,t) = G iiT (y,x,t) 
h 

+ 



+ 



4-7T 2 1 x — y I jo 

h 



dkk 



sin(/c|x — y|) sinh 2 ^/^) coth ^ i,k j 



\ ~ 2 

\ ^j,k^/ 5 k 



4vr 2 | 



\ dkki — 



y\ Jm f v ^,k^/ ; k 



sin(A;|x-y|)sin 2 (a; /ik t)coth(^^). (6.29) 



The first terms in Eqs. ( |6.28| ) and ( |6.29| ) are the two-point vacuum and thermal correlations 
( |6.23| ) and ( 6.24 ), respectively, before the quench. Therefore the remaining two terms de- 
scribe the effect of the quench. In particular, the second terms are dominant and rooted on 
the instability during the phase transition, which is missing in the field theoretical approach 
to equilibrium dynamics. Note that u; 2 k — u>j k 



mj + mj and tu 2 k + uj >k 



mf + m 2 , 

so the amplitudes of sin(/c|x — y|) decrease as 1/k 3 for very short wavelengths, hence short 
wavelengths contribute negligibly. However, there is a residual contribution from near the 
critical wavelength k c = rrif, which becomes much smaller than the second terms at later 
times and will not be considered any more. 

We wish to determine the size of domains from the second order phase transition of the 
instantaneous quench. At later times (m/t ^> 1) after the quench the dominant contribution 
to Eq. ( |6.29| ) comes from the second term, so one has approximately 



G 



h 



16n 2 r 



dk{ke 2Jj ^ t }sm(rk)F 1 (k), 



(6.30) 



where r = |x — y|, and 
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Fi(fc) 



'm 2 + m 2 f \ .ifPfoi*. 
1 coth 



(6.31) 



The function of k in the curly bracket of Eq. fl6.30| ) has a sharp peak at k$ = ^JmJ/2t, 
whereas Fi(k) is a slowly varying function. We employ the steepest descent method (see 
Appendix C) to obtain 



sin [ \ I m ^ v ) 
G fU;T (r,t) ~ G /c/ , T (0,t) V L^ ; ex P 



(6.32) 



2/ 



where 



(6.33) 



16?r 2 \ 2e^2t s _i V" u V 2t 

Therefore the size of domains grows according to the classical Calm-Allen scaling relation 



= J — ■ 

V m / 



(6.34) 



The scaling relation (|6.34j ) for the instantaneous quench confirms the result obtained in Refs. 



B. Finite Smooth Quench 

The instantaneous quench does not exhibit the essential spinodal behavior during the 
quench process. To see the dynamics of the second order phase transition one needs a finite 
quench period. Such a finite and smooth quench model is described by a field with the mass 
given by 

m 2 (t) = —ml — rngtanh^— (?71q > \m\\). (6.35) 

At earlier times t = — oo, the mass has the initial value 

m 2 = m\ = m\— m\ > 0, (6.36) 

and at later times t = oo, the final value 

m 2 = -m) = -(ml + m\) < 0. (6.37) 

Here r measures the quench rate: the large r-limit implies that the mass changes slowly 
from m 2 at t — — oo to — m 2 at t = +oo, whereas the small r-limit implies a rapid change 
of the mass. In particular, the (r = 0)-limit corresponds to the instantaneous change from 
m 2 to — m 2 at t — 0. That is, the instantaneous quench is a special case of the finite smooth 
quench model. To find the Fock space for each mode one needs to solve the classical equation 
of motion 
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ip k (t) + (k 2 - mj - m 2 tanh(-))^ k (i) = 0. (6.38) 

It should be noted that, as in the instantaneous quench model, long wavelength modes 
(k < k c = rrif) let the frequency change the sign at later times (t ^> r) 

wl(t) = k 2 - - ml < 0, (6.39) 

and suffer from the spinodal instability. Each long wavelength mode has a different quench 
time determined by u; k (£ k ) = 0. 

The solutions to Eq. ( |6.38|) are found separately for the stable modes and unstable 
modes. The stable modes (k > mj) have the solutions 



Mt) = C k e 2pk ' 2 Fi(/5+,k,/3_, k ;7k; -e 2t/T ), (6.40) 



where 



1 

T 



7 k = 1 - iru i;k , (6.41) 



with 



Wi jk = \fk 2 + m 2 , w/ jk = ^/A; 2 - mf. (6.42) 
Whereas the unstable modes (k < mj) have the solutions 

¥>k(*) = C k e 2 ^ 2 F 1 ( / 3 +jk , y 3_ jk ; 7k ; -e 2 ^), (6.43) 

where 

/3±,k = -^(^i,k ± £/,k), (6.44) 

with 



u> fM = ^/ m j - k 2 . (6.45) 

At earlier times (r <C —r) before the quench begins, both the solutions (|6.40|) and (|6.43|) 
have the same asymptotic form 

ViM = C k e- l ^\ (6.46) 
so the constant is fixed to satisfy Eq. 



c k = , 1 ■ (6.47) 
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1. During the Quench 



During the quench process (|t| < r), the asymptotic forms for the solutions ( |6.40[ ) and 
( |6.43| ) are unfortunately not available. Instead, one may expand the mass ( |6.35| ) to the 
linear order 



m 2 {t) 



t 

rrin[- 



(6.48) 



which is a good approximation as far as \t\ Cr. But we assume m§r 3> 1 so that the linear 
quench process is sufficiently long enough to allow the asymptotic analysis. Each mode of 
the scalar field then has the frequency squared 



1 2 2 

k — m 1 — m, 



•4) 



t k -t 



(6.49) 



where tk is the quench time for the corresponding unstable mode 



T 

ml 



(k 2 



m{) 



(6.50) 



Unless t| <C 1, the quench time occurs outside the valid regime for Eq. ( f[4D . 
So we restrict our attention to those unstable modes with <C r. The frequency ( |6.49|) 
corresponds to the special case of Eq. ( [4.17|) , in which li = l 2 = 0, to = tk, — U = T and 
io i = m . Then the solution in the intermediate regime before the quench, matching with 
the initial solution (|6.46[) , is given by Eq. (|4.18|) : 



¥W(t) = D 1 z 1 k /3 H[%(z k ) + D 2 zl /3 H[%(z k ), 



/3V 



where 



1/3 



(6.51) 



(6.52) 



The coefficients D\ and D 2 are given by Eq. (|4.20 ). After the quench time (t > tk) for each 
mode, the solution ( |6.51D is analytically continued for the unstable mode 



+ - 



D 



2 -1/3 



e^H^iz^+e^H^-iz^ 
■,,w - k ; ■ • ui/ 3 (— z5k) 



e*i "^^(^k) + e^tf 



(6.53) 



where 



t-t* 



1/3 



(6.54) 



The two-point thermal correlation function 
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G m ,T(r,t) 



2tt 2 Jo 



, 2 sin(/cr) 
dkk 2 



kr 



2tt 2 



,s m(kr) „ 



(*)Vm s ,k(t), 



(6.55) 



with r = |x — y|, is dominated by the unstable modes fl6.53|) , since the stable modes (|6.51|) 
oscillate rapidly and do contribute little. By using the asymptotic form for the solution 
( |6.53| ) in the regime r t ^> (r/mg) 1 ^ 3 



¥W,k(*) = 
one has approximately 




1 e Zk 



2tt ?y 6 

k 



(6.56) 



j_rj_^ v /3 s . 

8h m ^ Zk 



sin 



,2% 



(6.57) 



where Zi = 2moT/3. So the correlation function takes the form 

h 



G mu ,T{ r , t) 



16ir 2 r Jo 



dk{ke 2 "f< ki ) sin(kr)F n (k), 



where 



Fn(k) = — >' 



m y z k 



sin 2 (-m r + — ) cothf 



(6.58) 



(6.59) 



The function in the curly bracket is rapidly varying and has a peak at ko = {mo/^^/ri) 1 / 2 . 
Hence after applying the steepest decent method (see Appendix C), we finally obtain 



sin (&r) 

G mu>T (r,t) ~ G mUtT (0,t) ^° > 



cxp 



mo 



mo 



(6.60) 



where 



GW,t(0, t) 



ft 



7rm, 



3 \ 1/2 



\ l/Z 2 

\ (4 mo4+2 !^- r ) A /X / 



mo 



64n 2 m {(rt) 3 / 2 J 
We thus have shown the scaling relation for the domain size 



a/2 



•2Vrt 



(6.61) 



(6.62) 



The scaling relation (|6.62|) confirms, up to a numerical factor, the result in Ref. ||16|| . How- 
ever, the power-law is different from the Cahn- Allen scaling relation after the completion of 
quench. 
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2. After the Quench 

At later times (t ^> r) after the completion of quench, the solution ( J6.40 ) has the 
asymptotic form 



V/,k(*) 



2p k t 



2?UV. 



i,k 



r(7k)r(/3_, k -/3 + , k ) 
r(/3_ ik )r( 7k -/3 +)k ) 



e -2/3 + , ki /r 2j p i(/?+ik; 1 _ 7k + 1 _ ^ + fj +M . _ e -2t/r> 



r(7k)r(/3 +ik /3_ k ) 2/3 _ t/r „-2*/-n 

+ — vF7 « — V e 2-ri(.P-,k, l - 7k + P-k, J- - P+,k + P-,k, -e , 

1 (P+,k)l (7k - P-,kJ 



(6.63) 



The asymptotic form for Eq. (|6.43| ) is obtained by replacing p± )k by /3±, k . From the asymp- 
totic form of the hypergeometric function [ 42f , we find the asymptotic form for the stable 
modes 



<PfsM 



T(l - iuj iM r)T(-iuj f yr) 



and for the unstable modes 

1 



2hu itk r(l - i£(wi,k + w/, k ))r(-i|(a; ijk + w /)k )) 
1 r(l - zu;j )k r)r(ia; /ik r) 

2^j ik r(i - «IKk - w/,k))r(-z|(w i)k - ^ /ik )) 



r(i - i^, k ^)r(^/, k r) 

2^, k (-l)|(M,k - ^/, k )r 2 (-f (iw f)k - u) /)k )) 
T(l - za; i)k r)r(-a) /jk r) 



2hu iM (-i)l(^i,k + ^/, k )r 2 (-f (iw iik + & />k )) 



(6.64) 



(6.65) 



Now a few comments are in order. First, the coefficient of the positive frequency asymp- 
totic solution for each short wavelength mode 



2fuv 



(6.66) 



/,k 



leads to the rate for the initial vacuum to remain in the final vacuum 



|(0 k ,+oo|0 k ,-oo)p 



sufiV 



Y' (^ k + U SM 



sinh(7r™j jk ) sinh(7rra;/ ik ) 



(6.67) 



On the other hand, the coefficient of the negative frequency solution leads to the 

particle production rate ]35| 



1 - |(0 k ,+oo|0 k ,-oo>| ; 



sinh z 



sinh(7r™j jk ) sinh(7rrc<jy- jk ) 



(6.68) 



Second, when a)/ k r > 1, the coefficient of the decaying mode in Eq. ( |6.65|) can become 
infinite at 
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w/.kT = n, 



[n 



1,2,3,...) 



(6.69) 



because the gamma function has simple poles at these negative integers. For a rapid quench 
r < 1, there does not exist any integers that satisfy Eq. ( |6.69| ) . Hence this kind of resonance 
can happen only for a non-negligible r, i.e., for a very slow quench, which will be treated 
elsewhere. Contrary to the resonance at fl6.69D , the apparent singularity at £>/, k T = is 
removed by considering both terms in Eq. fl6.65| ): 



<Pfu*(t) 



r(l - iuiyr) 



2%u iM (-«|^i,k)r 2 (-i|o; iik ) 
1 T(l - iuivr) 



2?iu ik (-«l^i,k)r 2 (-i|o; iik ) 



r(-^ /ik r)e-^ k * + r(cD /ik r)e^ kt 
sinh(d) /)k t) 



^/,kr 



(6.70) 



After the completion of quench (t ^> r), the unstable modes ( |6.65| ) dominate the correla- 
tion function (|6.55|) over the stable modes (|6.64|) . In particular, the first term of Eq. (|6.65|) 
grows exponentially, so by using the asymptotic form in Appendix D, one has approximately 



8h 



/•k 



^i,k^/,k J sinh(7ru; iik ) [ 1 + ru)/ ik 



1 + |^/,k) + X^i.k 



n 2 



where £(n) is the Riemann zeta function and 



* = *-y(c(3)-i)(^ k + 



e V(C(2)-i)K k +^, k ) 

(6.71) 
(6.72) 



Note that + d)j k = mf + m^, so t lags by a constant in proportion to the cubic power 
of the quench period r. This time-lag is determined not only by the quench period but also 
by the initial and final coupling constants m-i and rrif. After applying the steepest decent 
method to the correlation function (see Appendix C), we finally obtain 



G fu , T (r,t)~G fU)T (0,t)- 



sin (\/? r ) 



m f , 

2t 



cxp 



rrifr 



8t 



(6.73) 



where 

Here 
F m (k) 



16tt 2 V 2e V 2t 



7T /m/N3/2 - , frnj 



(6.74) 



7TU Lk T 



^i,k^/,k / sinh(7ro;j )k ) 



1 +ra)/ |k 



1 2 



e V(c(2)-i)K k+ ^ k ) coth (^^). 

(6.75) 



Remarkably, the scaling relation of the domain size has the same form as Eq. ( |6.34| ) from 
the instantaneous quench 
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UAt) = .—. (6.76) 
However, there is a definite time-lag due to the finite quench as claimed in Ref. [fT3 j . The 



— 1 — — — o — i ii ii 

scaling relation ( |6.76| ) is robust, because only the asymptotic form of the exact solutions 
( |6.43| ) is used. 



VII. BACK-REACTION IN *7(l)-THEORY FOR PHASE TRANSITIONS 

The free scalar field model in Sec. VI does not describe a real system for the second 
order phase transition because the spinodal instability continues indefinitely. This model 
describes more appropriately an intermediate process of phase transition toward the spinodal 
line. In this section we consider |$| 4 -theory, in which there is a natural exit for the spinodal 
decomposition. The spinodal decomposition ends at the spinodal line, the regime in which 
the instability from the quench competes with the back-reaction from the |$| 4 -potential 
term. 

The ?7(l)-model for the second order phase transition is described by the Lagrangian 
density 

\ ( 12m 2 (t)\ 2 
C(t) = $*$ - V$* • V$ - - I <£>*$ + -Ai j , (7.1) 

where m 2 (t) is assumed to take either ( |6.21| ) or ( |6.35| ). Hence, the coupling parameter m 2 (t) 
of the quadratic term starts with the positive constant value mf at earlier times before the 
quench, keeps decreasing across zero at the moment of quench to, and finally reaches the 
negative constant value — mj at later times after the quench. So the $ = remains the true 
minimum of the potential until to- However, $ = is no longer the true minimum after t , 
but becomes a local maximum. The true minimum occurs at |$| 2 = 12m 2 (t)/A, and the 
system undergoes a second order phase transition. When the system is in the initial thermal 
equilibrium at earlier times, it is invariant under $ — > <3>e l61 and has a global L r (l)-symmetry. 
As the system undergoes the second order phase transition, the symmetry is broken and 
there occurs topological defects. 

^From the Lagrangian density (|7. 1|) follows the Hamiltonian density 



H(x, t) = n*n + V$* ■ V$* + m 2 (£)$*$ + ^($*$) 2 , (7.2) 

where we neglected a time-dependent c-number term. Since it has been shown in Sec. 
V that the coherent and coherent-thermal states in the LvN approach give the identical 
result as the mean-field and Hartree-Fock methods, we first divide the field into a classical 
background and quantum fluctuations and then quantize the latter according to the LvN 
approach. The field may be divided into the homogeneous classical background field and 
quantum fluctuations 

$(x,t) = c (t) + $ / (x,t) (7.3) 
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such that the classical background is a coherent state of the quantum field and the quantum 
fluctuations have symmetric states such as the vacuum, number and thermal states: 



($) = c (t), ($ / ) = 0. 



(7.4) 



The Hamiltonian density is then the sum of the classical background, quantum fluctuations 
and interactions 



W(x, t) = H c (t) +H f (x,t)+ W int (x, t) + 5H int {x, t), 



(7.5) 



where 



W 7 (x, t) = U* f U f + V$* f ■ V$/ + m 2 (t)®}$ f 



H mt (x,t) = -^($}$ / ), 

5H int (x,t) = vr c (n} +IL f ) + (V0 C ) ■ (V$} + V$/) + m 2 (t)0 c ($} + %) 

A f " ' ' ' * ^ 2</> c ($}$ / ) ($} + $/))• 



+^j^($; a + $j) + 2^($j + $ / ; 



(7.6) 



with 7T, 



C VC) 



As we are interested in the symmetric state of quantum fluctuations 



$/, so the expectation value of 5TLi nt vanishes and will not be considered any more, since 
only the terms involving ($^$j) do not vanish when one takes the expectation value with 
respect to the symmetric state. 

Now the field and $j and their conjugate momenta are decomposed into Fourier- 
modes according to Eqs. ( |6 4| ) and ( |6.5| ). The Fourier transform of the quartic term leads 
to 



/ rf 3 x($}$ s f = f[ / ^{<&M{fi«M6Qti + k 3 - k 2 - k 4 ). 



(7.7) 



A J (2vr) 3 ' 

The symmetric state further restricts the integral in Eq. ( |7. 7| ) to either (k x = k 2 , k 3 = k 4 ) 



or (k x = k 4 , k 2 = k 3 ), so Eq. (|7.7|) leads to 



( / d 3 x($}$ f ) 2 ) = 2 



_fk_ 

(2^)3 



2(<f>}!> / ) 2 . 



(7.8) 



Keeping the symmetric state in mind, we find the the Fourier transform of Tif + Timt 

d 3 k 



Ttf + 'Hint — J 



(2vr) 3 
A 



2 ■ 3! 



<(t)7T k (t) + (k 2 + m 2 (t) + -<P 2 c (t))K(t)Mt) 
d 3 k 



(2n) 



(7.9) 



Therefore the Hamiltonian for fluctuations is mode-decomposed into the sum of infinite 
number of coupled anharmonic oscillators. 
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We apply the method for quantum anharmonic oscillators in Sec. V to the Hamiltonian 
( [7.91) . According to the LvN approach, we first introduce two pairs of the annihilation and 
creation operators 



(7.10) 



and 



k h 



(7.11) 



We then require the creation and annihilation operators ( [7.10 ) to satisfy the LvN equation 
to obtain the equation of motion 



k 2 + m 2 (t) + ^(t) + 



A f rf 3 k! 



3! J (2vr) : 



<pu(t) = 0, 



(7.12) 



where the expectation value is taken with respect to the vacuum state, a symmetric state, 
as mentioned. The LvN equations for the operators ( |7.11|) are the complex conjugate of Eq. 

The Fock space for each mode can be found similarly according to Sec. V. We consider 
two symmetric states: the vacuum and thermal states. The vacuum state that is annihilated 
by all A k and 



A k (t)|0,t)=0, A k (t)|0,t> = 
leads to the expectation values 

d 3 k 



(7.13) 



($/$/)v 

(n* / n / )v 



(2vr) 3 
d 3 k 
(2vr) 3 



n 2 yi{t)y k {t) 



(7.14) 



The initial thermal state defined by the density operator for each mode 

' 1 



M*) = IIM*) = n )Y exp 



x^exp 



-(3hu i>k {A\{t)A k {t) + ^ 
(3%u> i!k (A£{t)Al{t) + i 



(7.15) 



leads to the expectation values 



Tr 



(L T }L T / ) T = Tr[p T (t)l>*<i> 



ft wmit) coth ( 2 , 



(2tt) 3 
d 3 k r 



h 2 tp h (t)<p k (t) coth( 



(7.16) 



37 



Then the effective Hamiltonian for the classical background with all the contributions from 
fluctuations in the initial thermal state is given by 



H c {t)=H c +{H int ) T 



vr 6 2 + 



m 2 (t) + 



A r d 3 k 



3\J (2tt) 



nY*(t)Mt) coth( 



(7.17) 



and that for fluctuations by 

H F = (Hf + Hint)? 

d 3 k 



(2tt) 3 
A 



A 



n 2 m)Mt) + (k 2 + m\t) + - } <p 2 c (t))h^(t)Mt) 



coth( 



2 ■ 3! 



^-^^ VkW^k(t) coth(— 



(7.18) 



The effective Hamiltonian from the initial vacuum state is the zero-temperature limit (/5 — ► 
oo), i.e., coth(/3^o;j i k/2) — > 1. 

On the other hand, in the effective action method of Sec. V.B, the tp^ is a parameter 
that will be determined by the Hamilton equations. By writing ip^ in the polar form 



Ck(t) __ 



(7.19) 



and by introducing = £k, the effective Hamiltonian for the k-mode ( |7.18| ) can be rewritten 
as 



#T,kOO = ^coth( 



where 



ll4(t ) = k . + m . (t) + ^| W + »/^ 



coth( — — Jf/v,k(^J- 



(7.20) 



Ck(*) coth( 



(7.21) 



The Hamilton equations are 



9 



#T,k(*) 



^Ck 

^PCk L ^coth( 



(9 



#T,k(*) 



<9Ck L frcoth( 



-Wk(*)Ck + 



4C; 



(7.22) 



These Hamilton equations are identical to Eq. ( 7.12Q , because the effective action method 
is equivalent to the LvN approach as shown in Sec. V.B. The Hamilton equations of motion 
for the classical background are given by 



m 2 {t) + 4 



' /: '' k ^VkW^k(t)coth(^; 



(2tt 



<t>c{t) 



A 



2 ■ 3! 



0. 



(7.23) 
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The different quantum states given in Sec. V can be chosen to describe the different 
processes for the phase transition. We assume that the system starts from either the ther- 
mal or coherent-thermal states before the quench and evolves according to the functional 
Schrodinger equation during and after the quench. In this sense the process is completely 
determined once the initial quantum state is prescribed. In the first case of the symmetric 
thermal state, each mode has the zero expectation value for the position and momentum, 
but its dynamics is still governed by Eq. ( [7.12[ ) with C = 0. In the second case of the 
coherent-thermal state, the classical field C , which is a coherent state of the homogeneous 
field, plays the role of an order parameter and is influenced by thermal fluctuations of 0k. 
As the first case is a limit of the second one, we shall first focus on the coherent-thermal 
and then obtain the thermal state result by taking the limit 4> c (t) = <j> c (t) = in the end. 

It is very difficult to solve analytically the equations of motion ( |7.12 ) and ( |7.23| ). At best 
we have to rely on the adiabatic solutions that can be found in some important physical 
regimes. We assume that far before the quench (t — > — oo) the system starts from a thermal 
fluctuations around the order parameter with 



-oo 



0. 



-oo 



0. 



(7.24) 



As C is a classical field, the uncertainty principle does not prohibit us from even taking the 
limit C (— oo) = C (— oo) = 0, which leads to the symmetric thermal state. The initial data 
for C take very small values with high probability by some distribution function. For the 
sake of simplicity we consider the instantaneous quench first. Before the quench time, m 2 (t) 
takes the initial constant value mf and C remains close to zero, so changes very little 



during the evolution. The solutions to Eq. ( 7.12 ) are approximately given by 



2Mh 



and Eq. (|7.21|) yields the gap equation 



mf + k 2 + — 



(7.25) 



i.k 



Xh r d 3 k l 



3! J (27r) 3 2fi ijk 



coth( 



(7.26) 



The infinite quantity that appears in Eq. (|7.26j ) will be absorbed by the bare coupling 
parameters m 2 and A to result in the renormalized ones and a finite equation for f2? k [H. T2 
The solutions ( [7.251 ) hold till the quench time t = 0. 



But after the quench time, m 2 (t) changes to 



-m 



/' 



and therefore, the long wavelength 



modes become unstable and the spinodal instability begins. The long wavelength solutions 
are given by 



l-i ^^ sinh( / n f Jt)) +cosh( / fi fk (t)) 



(7.27) 



and the short wavelength solutions by 



V/s,k(*) 



2TiVL 



—i 



i.k 



i,k 



%k(t) sin( io^ k(t))+cos( yo^ k(t)) 



(7.28) 
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where 



m 



f 



k2 ~ ^ 



Xh r k ^ d 3 k coth( ., 



Aft 



(2tt) 3 2n iik 
rf 3 k coth(^i 



3! Jk A (2vr)3 2fi iik 




?£)) 2 + 1 } ^(/;%, k W) + i 



fy,k(t) 



l| sin(^ % k (t)) + l , (7.29) 



and fiy k (t) = — f2j k (t). In contrast with the free scalar model in Sec. VI, not only the 
duration of the instability but also the band of unstable modes decrease in time due to the 
back-reaction ft 2 <^£(£)<^ k (t) of the self-interaction and 4> 2 c (t) of the classical background. The 
classical background obeying Eq. Q7.23 ), though it remained around the initial vacuum under 
thermal fluctuations before the quench time, begins to roll down from the false vacuum to 
the true vacuum at first largely due to the exponentially growing unstable modes and then 
to the self-interaction. The competition between — mj and ($^$^) together <fi 2 c determines 
the spinodal line, beyond which the instability stops and the fluctuations begin to oscillate 
around the true vacuum. Therefore the self-interacting phase transition model has a natural 
exit to the spinodal instability [2^ . 

We now turn to the case of symmetric thermal state. Once the initial condition is 
prescribed such that C (— oo) = C (— oo) = 0, the classical field 4> c (t) remains in the false 
vacuum even during the quench and has the trivial solution <p c (t) = for all times. Hence 
Eq. (|7. 23|) is identically satisfied, and the dynamics of the second order phase transition 



is entirely governed by Eq. ( 7.12|) for quantum fluctuations, which extends the free scalar 
field model considered on Sec. VI. Though each mode of quantum fluctuations has the 
zero expectation value, the Wigner function becomes sharply peaked around its classical 
trajectory as the quench proceeds Even without the classical field C the quantum 

contribution from the self-interaction in Eq. (|7.12|) still prevents the unstable modes from 
growing indefinitely and provides a natural exit to the spinodal instability. This implies 
that the quantum dynamics of the second order phase transition is classically correlated 
and exhibits most of the essential points described by the order parameter under thermal 
fluctuations. 

Finally we comment on the formation process of topological defects. The topological 
defects formed from the second order phase transition in this paper are domain walls. The 
correlation of domain walls can be determined by the two-point thermal correlation function 



G T (y,x,t) = <6*(y,t)6(x,t)) T 



d 3 k r 



(27T) 



h 2 ip* k {t)ip u {t) coth( 



Q ik-(x-y) 



(7.30) 



where a>j >k = cj k (t = — oo). After the quench, the two-point correlation function is dominated 
by the unstable modes (pf v . The scaling behavior of two-point correlation function in Sec. VI 
holds still before reaching the spinodal line. That is, the size of domains grows according to 
the power law t 1//4 during the quench and the Cahn-Allen relation t 1//2 after the completion 
of quench. However, the domains can not grow indefinitely due to the back-reaction of 
the self-interaction and the classical background. As each unstable mode reaches further 
the spinodal line, its solution stops exponential growing and oscillates. The behavior of 
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the two-point thermal correlation function changes from that in Sec. VI, and one would 
expect a different scaling relation for the domain size, probably with a small power than the 
Cahn-Allen relation. 

VIII. CONCLUSION 

In this paper we have elaborated the recently introduced Liouville-von Neumann (LvN) 
approach to describe properly the time-dependent nonequilibrium systems. The systems 
interacting directly with environments or undergoing phase transitions are such nonequilib- 
rium systems. These systems are characterized by time-dependent coupling parameters and 
their true nonequilibrium evolution deviates significantly from the equilibrium one when 
their coupling parameters differ greatly from their initial values. In this case the systems 
evolve completely out of equilibrium. For that purpose there have been developed many 
different methods such as the closed time-path integral method, sometimes in conjunction 
with the large iV-expansion, mean-field, Hartree-Fock method. The LvN approach devel- 
oped in this paper is a canonical method that unifies the functional Schrodinger equation 
for the quantum evolution of pure states and the LvN equation for the quantum description 
of mixed states of either equilibrium or nonequilibrium. Because the LvN approach shares 
all the useful techniques with quantum mechanics and quantum many-particle systems, it 
turns out to be a powerful method for describing time-dependent harmonic oscillators and 
anharmonic oscillators, and provides a rigorous and systematic method for describing time- 
dependent phase transitions. 

By applying the LvN approach to time-dependent harmonic oscillators, we have found 
exactly the nonequilibrium quantum evolution evolving from various initial states such as the 
vacuum, number, coherent and thermal states. In this case the LvN approach is based on two 
operators, the so-called the annihilation and creation operators, that satisfy the quantum 
LvN equation, so it is straightforward to construct the Fock space of number states and 
the density operator according to the standard technique of quantum mechanics. We have 
thus obtained the density operator in terms of the classical solution, and by using the exact 
wave functions for number states, have been able to find the explicit form of the density 
matrix. In particular, the density matrix provides us with a criterion on nonequilibrium 
vs. equilibrium evolution. Moreover, the LvN approach has been applied to the time- 
dependent inverted harmonic oscillators, which can be regarded as models for second order 
phase transitions. For time-dependent anharmonic oscillators, we have found approximately 
the nonequilibrium evolution of the symmetric Gaussian, coherent and thermal states at the 
lowest order of the coupling constant of the quartic term. It has been shown that the LvN 
approach is equivalent to the effective action method and to the mean field or Hartree-Fock 
method. 

Finally we have applied the LvN approach to the systems undergoing the symmetry 
breaking second order phase transition. In particular, due to the quench the coupling pa- 
rameters change the sign during the evolution. As field models we have studied a free massive 
scalar field with an instantaneous and a finite smooth quench. By applying the LvN ap- 
proach to this symmetry breaking system we have found the two-point vacuum and thermal 
correlation functions. It has proved that the spinodal instability leads to the t 1 / 4 -scaling 
relation for domain sizes during the quench and the classical Cahn-Allen relation after the 



41 



completion of quench. The Calm-Allen scaling relation confirms the result for the instan- 
taneous quench model in Refs. 0,0. One prominent feature of the finite smooth quench 
model is the time-lag occurring at the cubic power of the quench period in the Cahn- Allen 
scaling relation after the completion of quench. The inclusion of a self-interacting term shuts 
off the spinodal instability after crossing the spinodal line and gives rise to a natural exit 
for the spinodal decomposition. Not treated in detail in this paper is the very slow quench 
effect, which may show a transient resonance of decaying solution of long wavelength modes 
and will be addressed in a future research. 
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APPENDIX A: HARMONIC OSCILLATOR WAVE FUNCTIONS 

The vacuum state of the time-dependent oscillator is annihilated by a(t): 



In the coordinate representation 



a(t)\0,t) = 0. 



* (q,t) = (q\Q,t), 



Eq. (|A1|) becomes 



u — — mu q 

i aq 



*o(?,*) = 0. 



We thus obtain the normalized wave function for the vacuum state 

1/4 

exp 



i 



2nh 2 u*u 



i mu 2 



(Al) 



(A2) 



(A3) 



(A4) 



The wave function for the nth number state is obtained by applying the creation operator 
a'(t) n-times 



V n (q,t) = -±=(a^t)) n y (q,t). 
By making use of the relation 

(a f (t)J *(g,t) = (-ftuj ez-- q (—) ( e~^~^ q t) 
and the definition of the Hermite polynomial 

H n (x) = 



(A5) 



(A6) 



(A7) 



42 



we obtain the wave function 



1/4 



2irU 2 u*u) V2 n n! \u 



H n (x)exp 



i mu 2 
2~hv* q 



where 



x 



q 



V2h 2 u*u 

Here we have also used the wronskian (13.51). 



APPENDIX B: DENSITY MATRIX 



In the coordinate representation, the density operator defined by 



- / N 1 

PT[t) = 

Zj N 



-/3hu (N(t)+±) 



where is the partition function, becomes 



By substituting (|A8[) into (|B2|), we obtain 



If 1 



1/2 



Z N \27rh 2 u*uJ ^ 2 n n\ 



£ —M n {x')H n {x)e 



ftp,, . ( rrt _i_ M * m £ m U„£ 



Following Kubo's method | 40fl , we rewrite the product of Hermite polynomials as 

\ n — zf — 2ix' z\ — z^ — 2ixz2 



H n {x')H n {x) = -e x " + '~ 



and sum over n to obtain 



oo roc 



dz l dz 2 (2iz l ) n (2iz 2 ) n e~ 



1 



Zj\f \2irfi 2 u*u 



1/2 _/3R^0 

e 2 



7T 



exp 



-'2 



U 



{x 2 + x 2 ) + - — < —q - - -q~ 
2 n \ u u 



x 



dz\dz 2 exp 



oo J — oo 



— — 2ix'z\ — z\ — 2ixz 2 — 2z\z 2 e 



After doing the integral and using the identity 



2 i m fU u 
X ~ 9¥W ~^*' lJ 



we finally obtain 



Airh u*u sinh(/%u;o) 



1/2 



exp 



—Hm—\n(u*u)(x' 2 — x 2 ) 



x exp 



\\{x' + x) 2 tanh(^) + (x> -*) 2 coth(^p)' 
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In the special case of the time-independent oscillator we recover the classical result by Kubo 
4Cj by substituting the complex solution 



u{t) 



-iujQt 



into Eq. 



APPENDIX C: STEEPEST DECENT METHOD 



The integral appearing in the two-point correlation function has the form 

1=1 dxxe'" 1 * 2 sm(yx)F(x) (CI) 

where F(x) is a slowly varying function. Note that xe~ lx2 is a highly peaked function that 
varies rapidly. We let 



where 



e 9{x) = xe~ lx2 



g(x) = — 7X 2 + ln(x). 



(C2) 



(C3) 



We expand g(x) in a Taylor series and truncate it up to the quadratic term around the 
maximum point Xq = -7= 

g(x) « g(x ) - 2 7 (x - x ) 2 , (C4) 

and rewrite the integrand as 

xe^ sm(yx)F(x) « x e-^F{x )e~ 2 ^ x - v ' ] 
= xoe~ lx2 °F(xo) \ exp 



\2 e 



2i 



of ■ y \* 



O iyxo 



2i 



+ exp 



-27 ( 



x — xq + i 



y_ 
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(-l)e- tyxo ) _ vl 

>e 87. 

2i I 



(C5) 



The Gaussian integrals contribute equally, so we obtain 

1/2 



7T 



(C6) 
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APPENDIX D: ASYMPTOTIC FORM FOR UNSTABLE MODES IN A FINITE 

QUENCH 



To find the contribution to the two-point functions from the unstable growing 
for a finite quench (u)/, k T < 1), we need to evaluate r(u;/ i k) and r(|(u)/ )k — iuji^)). 
gamma functions are rewritten as 



modes 
These 



r(i + a /lk T) 

r(i + f(d)/,k -iuj,k)) 



r(-(%,k-^ ik )) 



We further make use of the expansion formula |42| 

oo ,, 

lnr(l + z) = - ln(l + z) + (1 - 7)2 + EM)" [C(n) - l] -, (\z\ < 2), 

where 



ra=2 



7 = lim 



Ei- ln ( m ) 



./c=i 



0.5772156649 



is the Euler's constant and ((n) is the Riemann Zeta function 

00 1 



Here 



k=l 



T 



We now expand the gamma functions up to the cubic power of r: 



r(l + a) /)k r) 



1 



1 + £>/,kT 
1 



exp 



;i-7)(w/ > kT) + i[C(2)-l]fer) 



[C(3)-l](^ )k rr + 0(r 4 ) 



and 



T 



r(i + -(%, k -^i,k)) 



|1 + |(o)/,k - ^i,k) 



exp 



;i-7)(a; /)k r) 



+ ^[C(2) - l](^) 2 {(^/, k - ^, k ) 2 + (a} /)k + ^, k ) 2 } 

" i[C(3) " l](~) 3 {(%,k - ^,k) 3 + - ^, k ) 3 } + 0(r 4 



(Dl) 



(D2) 



(D3) 



(D4) 



(D5) 



(D6) 



• (D7) 



Therefore it follows that 
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r(l+a> /lk r) 



P(l + |(o)/,k - ioJi,-k)) 
x exp 



(» 



2 r 2 2 



1 + &f,VT 

Si 



i[C(2) - l]r 2 (^ k + wj k ) - i[C(3) - l]r 3 a) /)k (4 k + + 0(r 4 



We finally get 



r (^/,kT) 


2 


"Tfe + ^k)" 


2 


r(i + i^/,k) 2 + ^ki 2 


r 2 (f (o)/, k - ^i,k)) 











x exp 



i[C(2) - l]r 2 (^ k + c 2 k ) - i[C(3) - l]r 3 o) /)k (4 k + c 2 k ) + G(r 4 ) 



The last relation necessary for the two-point correlation function is fl42 

2 



r(l - za» i)k r] 



7r^ ik r 



sinh(7ra;j )k r) 



(D8) 



(D9) 



(D10) 
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